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Abstract. We show how the argument exploited by Galaverni & Sigl in 
Phys. Rev. Lett., 100, 021102 (2008) (see also arXiv:0708.1737) [1] to constrain Lorentz 
invariance violation (LV) using Ultra-High-Energy photon non observation by the 
AUGER experiment, can be extended to QED with Planck-suppressed LV (at order 
0{E/M) and 0{E'^/M'^)). While the original constraints given by Galaverni & Sigl [1] 
happen to be weakened, we show that, when used together with other EFT reactions 
and the expected detection of photons a± E > 10^^ eV by AUGER, this method has 
the potentiality not only to basically rule out order 0{E/M) corrections but also to 
strongly constrain, for the first time, the CPT-even 0{E^/M'^) LV QED. 



1. Introduction 

The interest in possible high energy violations of local Lorentz Invariance (LI) has grown 
in recent years, encouraged by the flourishing of observational tests and the availability 
of experimental data. Theoretically, hints of Lorentz Violation (LV) arose from various 
approaches to Quantum Gravity (QG) [SlEliaElEllTlElE]. 

While it may seem hopeless to search directly for effects suppressed by the Planck 
scale, it is possible that even tiny corrections can be magnified to measurable ones 
when dealing with high energies, long distances of propagation or peculiar reactions 
(see e.g. p^Of [TT]). However, in order to do so, a complete theoretical framework, within 
which it is possible to calculate reaction rates and to describe the particle dynamics, is 
needed. 

This is the case of Effective Field Theory (EFT) with LV operators. We shall deal 
here with modified QED via non-renormalizable, Planck suppressed LV operators (the 
analogue theory with renormalizable operators being already severely constrained |10j). 
It has been shown fl2[ [13] that the addition of the two lowest order non-renormalizable 
LV operators (mass dimension 5 and 6 respectively) to the effective Lagrangian of QED 
leads to the following high-energy modified dispersion relations (MDR) (since the LV 
correction is proportional to p" we call them n— MDR with n = 3 for dimension 5 
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operators and n = 4 for dimension 6 ones) 



(1) 
(2) 



where ([T]) refers to photons while ([2]) refers to fermionslj. (In the following we will assume 
M to be comparable to the Planck mass Mpi ~ 1.22 x 10^^ GeV.) The constants 
and r^j!^'' indicate the strength of the LV and take values on the whole real axis. In ([1]) 
the + and — signs denote right and left circular polarization, while in ([2]) they indicate 
opposite helicity fermion states. 

A crucial difference between the n = 3 and n = 4 cases is the fact that the former 
is characterized by LV terms which break CPT invariance, while the relevant ones for 
the 4— MDR are CPT even [13j. This difference implies that for the 3— MDR there is an 
effective breaking of the symmetry between the two helicity states of the photon. Indeed 
one finds that if = -^i^^ = i'^^\ while if = if = i^^^ |I3]. On the other hand, it can 
be shown that the coefficients of electrons and positrons are related as r/^. = {—)"-ri'f, 
exploiting the argument given in [13] . 

Since suitable powers of the suppressing scale Mpi have been already factored out 
in eq.([Il E]), natural values of the LV, dimensionless coefficients in ([1], [2]) are expected 
to be 0(1). 

From a purely logical point of view, it could seem unreasonable to study 0{E'^ / M"^) 
LV corrections (those leading to 4— MDR), as these will be always subdominant with 
respect to those at 0{E/M) (those leading to 3— MDR). However, the reasons for this 
interest are both empirical and theoretical. 

On the observational side, the LV parameters rj}^ are presently constrained to be 
less than O(10~^) at 95% confidence level (CL) by a detailed analysis of the synchrotron 
component of the Crab Nebula broadband spectrum [15], while the constraint |^*^'^^| < 
10"'' is obtained by considering the absence of vacuum birefringence effects in the 
propagation of optical/UV polarized light from Gamma- Ray Bursts [T6] . 

On the theoretical side, a reasonably good motivation for focusing on 0(i?^/M^) 
LV corrections is related to the so called "naturalness problem" [13]. In fact, it is generic 
that, even starting with an EFT with only mass dimension 5 or 6 LV operators for free 
particles, radiative corrections due to particle interactions will generate lower dimension 
LV terms without introducing any further suppression [T7|. Hence extra LV terms in 
p and p"^ will be generically dominant on higher order LV ones and lead to extremely 
stringent constraints on the dimensionless LV coefficients i^"'\r]^"-\ 

However, it has been shown [18] that if the theory includes SuperSymmetry 
(SUSY), then dimension 3 and 4, renormalizable, LV operators are forbidden. As a 
consequence, renormalization group equations for Supersymmetric QED with dimension 
5 LV operators a la Myers & Pospelov were shown not to generate lower dimensional 
operators, if SUSY is unbroken. 

I Actually, 4— MDR contains an extra term, proportional to p^, due to CPT even dimension 5 
operators [l3]. However this term is suppressed by me/A/pi. 
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SUSY soft-breaking does lead to LV terms in p and p^, however characterized 
only by a suppression of order m^/Mpi (n = 3) or (ms/Mpi)^ {n = 4), where 
~ 1 TeV is the scale of SUSY soft breaking [ISj. Nonetheless, given the present 
constraints, dimension 5 LV operators would induce dimension 3 ones which are 
already tremendously constrained. Hence, in the n = 3 case one would have to 
require unnaturally small LV coefficients in order to have a viable model. On the 
contrary, if n = 4 then the induced dimension 4 terms are suppressed enough to be 
compatible with current constraints without requiring , i]^^^ much less than one, 
provided nis < 100 TeV. Therefore, missing an alternative "custodial symmetry" for 
LV with respect to SUSY, QED dimension 5 LV operators seem problematic, while 
dimension 6 LV, CPT even, ones are favored. 

At present, a clear general argument, as to why dimension 5 LV operators should 
not appear, is missing. However, if we assume, together with SUSY symmetry with 
m<j < 100 TeV, also CPT invariance for the Planck scale theory, then not only dimension 
3 and 4, but also dimension 5, CPT odd, LV operators would be forbidden and only 
CPT even, dimension 6 ones would appear in the effective Lagrangian [13] 

Given this state of affairs, we shall provide here constraints on both dimension 5 
and dimension 6 operators separately. Of course, in order to do so for dimension 6 
operators one would need to explore much higher energies that those associated to the 
above cited constraints on n = 3 dispersion relations. This is why we shall focus here 
on the physics of ultra-high-energy cosmic rays (UHECR). 

It has been recently pointed out [Ij that if Lorentz symmetry was violated, then the 
absorption of ultra-high-energy (UHE) photons {E > 10^^ eV) on the Cosmic Microwave 
Background (CMB) and the Universal Radio Background (URB) could be forbidden, 
thus leading to large photon fluxes reaching Earth. This would violate limits put by 
current experiments [191 120] on the photon fraction in UHECR. Hence, very strong 
constraints |^(^)| < 10"^^ and ^(^^ > -10"^ were claimed [1] . An underlying assumption 
in [1] is that t]*^"-* ~ in order to prevent competing reactions with respect to photon 
absorption. This is an important limitation, as we will show below. 

The main result of this work is to extend the idea given in [Ij to the full LV QED 
framework described by equations ([1]) and ([2]). While the original constraints [1] are 
weakened, we shall see that, when used together with other EFT reactions, this method 
has the potentiality not only to basically rule out the n = 3 case but also to strongly 
constrain, for the first time, the CPT-even (hence possibly theoretically favored) n = 4 
LV QED. 

This paper is structured as follows: in section [2] we briefly describe the origin of the 
Ultra- High- Energy photons and explain the argument by flj. In section [3] we describe 
the generalization of the previous argument in Effective Field Theory LV. In section H] 
we discuss our results and draw conclusions. 

§ It is however important to stress that the SUSY LV operators considered in [18] do not lead to 
dispersion relations of the form presented here. 



GZK photon constraints on Planck scale Lorentz violation in QED 



4 



2. UHE photons and LV 

UHE photons originate in the interactions of UHECRs with the CMB, leading to 
the production of neutral pions which subsequently decay into photon pairs. Pion 
production occurs only if the interacting UHECR energy is above i?th — 5 x 
10^^ (ciJh/l.S meV)^^ eV (cjf, is the target photon energy). Hence, it has long been 
thought to be responsible for a cut off in the UHECR spectrum, the Greisen-Zatsepin- 
Kuzmin (GZK) cut-off |2l]. 

Experimentally, the presence of a suppression of the UHECR flux has been 
confirmed only recently with the observations by the HiReS detector [22] and the Pierre 
Auger Observatory (PAO) [23]. Although the cut off could be also due to the finite 
acceleration power of the UHECR sources, the fact that it occurs at the expected 
energy favors the GZK explanation. The results shown in [21] further strengthen this 
hypothesis. 

The PAO and the Yakutsk and AGASA experiments also imposed limits on the 
presence of photons in the UHECR spectrum. In particular, the photon fraction is less 
than 2.0%, 5.1%, 31% and 36% (95% C.L) aX E = 10, 20, 40, 100 EeV respectively 
[T9l [20] . Although its theoretical computation is quite uncertain and depends on many 
unknowns related to source and propagation effects [19], it is established that photons 
are mainly attenuated by pair production onto CMB and URB. 

However, pair production is strongly affected by LV. In particular, the (lower) 
threshold energy can be slightly shifted and in general an upper threshold (a finite energy 
above which pair production is no more allowed by energy-momentum conservation) can 
be introduced [26]. Therefore, if the upper threshold energy happens to be lower than 
10^^ eV, then UHE photons are no more attenuated by the CMB and can reach the 
Earth constituting a significant fraction of the total UHECR fiux, thereby violating 
present experimental limit^ [Ij. 

However, this argument is not stringent enough to cast constraints on LV in EFT, 
because in this framework two competitive processes, forbidden in LI physics, are allowed 
and can effectively dump the photon flux: photon decay in vacuum and photon splitting 
(7 — i> A^7). In [1] the special case r^^") ~ and < was considered, in order to 
prevent these extra processes. In the following we will study the full parameter space. 

3. LV reactions 

In order to perform a consistent analysis in LV EFT, we have to consider three processes 
related to photon attenuation: pair production, 7-decay and photon splitting. 

Pair production This well known process occurs whenever the center-of-mass energy 

II This conclusion could be evaded if the GZK process was not effective. However, the large mass 
difference between pions and electrons implies that, at comparable energies and LV coefficients, the 
GZK reaction must be much less affected than pair production. Moreover, it can be shown [1] that LV 
does not affect the kinematics of 7r° decay. 



GZK photon constraints on Planck scale Lorentz violation in QED 



5 



of the 77 system is sufficient to produce a e^/e" pair, i.e. it is larger than 2me. 
This condition corresponds to A; > fcth = fnl/^^b- 

In the following, we exploit the above mentioned relation r}^ = {—)""ri'^ , as on 
average the initial state is unpolarized (see however [25] for further discussion on this 
point). This is justified as 7r° decay produces photons with opposite polarization 
and the CMB is unpolarized on average. Moreover, the interaction at threshold 
must occur in 5- wave, because the particles' momenta need to be aligned [26] . 
Nevertheless, it is possible that, if the S'-wave channel is forbidden, higher partial 
mode interactions occur in off-threshold configuration. In this case, however, the 
reaction rate is suppressed by partial mode suppression and because only suitably 
polarized initial states can contribute to it. 

Within this framework, and exploiting energy-momentum conservation, the 
kinematics equation governing pair production is the following [Hj 



where ^ = ^W/M""^ and fj = r/("')/M"'~^ are respectively the photon's and 
electron's LV coefficients divided by powers of M, < ?/ < 1 is the fraction 
of momentum carried by either the electron or the positron with respect to the 
momentum k of the incoming high-energy photon and uJb is the energy of the target 
photon (we will assume in the following Ub = ^^cmb — 6 x 10^^ eV and will not 
consider pair production onto the URB, as our main conclusions can be drawn 
using just CMB). Note that the symmetry of ([3]) under the exchange ^ e~ is 
manifested in its symmetry under y ^ 1 — y and fj ^ {~)"'fj. 
Pair production can be severely affected by LV. In particular, it has been shown 
|26j that, rather surprisingly, not only the threshold energy fcth is modified, but 
also an upper threshold is introduced. Physically, this means that at sufficiently 
high momentum the photon does not carry enough energy to create a pair and 
simultaneously conserve energy and momentum. However, an upper threshold can 
only be found in regions of the parameter space in which the 7-decay is forbidden, 
because if a single photon is able to create a pair, then a fortiori two interacting 
photons will do [26] . 

The structure of the lower and upper thresholds for n = 3,4 has been studied 
in [26] if 77+ = ?7_. The same kind of analysis can be extended to the full EFT 
case. However, being the computation rather cumbersome, we shall evaluate it 
numerically. The structure of the constraint is different depending on n. If n = 3, 
since = —C- and the constraint is imposed on both left and right polarized 
photons, it is symmetric with respect to ^ —C,- If n = 4 such a symmetric 
structure is lost. 

7-decay While forbidden in LI theory, this reaction is allowed in a LV framework if 
the photon energy is above a certain threshold. The latter can be easily derived by 
solving the relative energy-momentum conservation equation which can be readily 




(3) 
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inferred from ^ by noticing that it corresponds to the hmit cjf, — ^ 0. 

Photon splitting This is forbidden for ^^"^ < while it is always allowed if ^'•"^ > 
|26j . When allowed, the relevance of this process is simply related to its rate. The 
most relevant cases are 7 — > 77 and 7 87, because processes with more photons 
in the final state are suppressed by more powers of the fine structure constant. 
The 7 — > 77 process is forbidden in QED because of kinematics and C-parity 
conservation. In LV EFT neither condition holds. However, we can argue that this 
process is suppressed by an additional power of the Planck mass, with respect to 
7 87. In fact, in LI QED the matrix element is zero due to the exact cancellation 
of fermionic and anti-fermionic loops. In LV EFT this cancellation is not exact and 
the matrix element is expected to be proportional to at least (^-E'/Mpi)^, p > 0, as 
it is induced by LV and must vanish in the limit Mpi — > 00. 

Therefore, we have to deal only with 7 — 87. This process has been studied in 
[26l [27] . In particular, in [27] it was found that, if the "effective photon mass" 
= ^i?"/Mpf^ <^ rrig, then the splitting lifetime of a photon is approximately 
^n=3 ^ 0.025^^^/^^ (50 leY / E^)^^ s, where / is a phase space factor of order 1. 
This rate was rather higher than the one obtained via dimensional analysis in [26] 
because, due to integration of loop factors, additional dimensionless contributions 
proportional to enhance the splitting rate at low energy. 

This analysis, however, does not apply to our case, because for photons around 
10^^ eV ^ ml if ^^^^ > 10"^'' and ^"^^^ > 10"^. Hence the above mentioned 
loop contributions are at most logarithmic, as the momentum circulating in the 
fermionic loop is much larger than mg. Moreover, in this regime the splitting rate 
depends only on m^, the only energy scale present in the problem. 
We then expect the analysis proposed in [26] to be correct and we infer that the 
splitting time scale at ~ 10^^ eV is larger than the propagation one (100 Myr 
for GZK photons) if ^^^^ < 0.08, while if n = 4 it is well above 100 Myr even for 
e(^) ~ 0(1). 

4. Results and conclusions 

We have demonstrated so far that only 7-decay and pair production are relevant to our 
analysis. By considering these processes three kinds of constraints are possible. 

On the one hand, since at present we have stringent upper bounds on photon fluxes 
up to 10^° eV [20] , then any upper threshold energy introduced by LV in pair production 
must be larger than this figure. This leads to the constraint represented by the black 
thick solid lines in Fig. [H where the allowed region is obviously the one including the 
origin and in the case n = 8 is the intersection of the upper threshold allowed region 
with the ones allowed by already existing constraints (red lines). We confirm the claim 
by pLj that, for r/(3,4) ^ g, < 10"^^ and > -10"^ However, Fig. □ shows that 
this is a rather special (and favorable) case. 

On the other hand, if some photons were detected above 10^^ eV, then it could be 




Figure 1. Left panel: n = 3 LV. Right panel: n — 4 LV. Constraints from the absence 
of pair production upper threshold. The best constraints to date are shown in red, if 
they exist. The allowed region includes the origin and corresponds to the intersection 
of the regions bounded by the red and black lines. 
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Figure 2. Left panel: n = 3 LV. Right panel: n = 4 LV. 7-decay threshold structure. 
The best constraints to date are shown in red, if they exist. 



deduced that the threshold energy for 7-decay is larger than this energy (once allowed, 
photon decay is basically instantaneous [H]). Indeed, the PAO will reach the required 
sensitivity to probe such theoretically expected fluxes [28] within the next few years 
[19]. In this case the allowed region for n = 3 would be the "clepsydra" delimited by 
the black solid line in Fig. [2] and by the two horizontal red lines corresponding to the 
birefringence constraint. As it can be inferred from Fig. [21 for n = 3 the combination 
of the above mentioned constraints would cast the very strong bound \ri^^''\ < 10~^. 
Conversely, no significant limit would be placed for n = 4. 

Finally, we notice that these two methods do not in principle exclude each other 
(although the physical effects are mutually exclusive). This is the case if a lower limit to 
the photon flux at 10^^ eV is imposed by the experiments (thus implying the absence of 
7-decay), while the upper limit at 10^° eV is confirmed. The constraint obtained in this 
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case would be very strong, as the allowed region is given by the intersection of the two 
regions bounded by the green dot-dashed and the black solid lines in Fig. [31 Figure [3] 



-6 -8 -10 -12 -14 -If -1 



16 -14 -12 -10 -8 -6 



-4 -2 



Crab 



2 









-2 




-4 




-6 




2 0-2-4 -6 


-8 
-10 

-8 -10 


-10 -8 -e -4-2 2 




-10 
-8 




y/^ -6 




y< -4 




-2 








--''^ 2 





Figure 3. Left panel: n = 3 LV. Right panel: n = 4 LV. The LV parameter space is 
shown. The current best constraints (when they exist) are drawn in red. Black solid 
lines represent values of (77, ^) for which the 7-decay threshold k^-dec — 10^® eV. Dot- 
dashed, green lines indicate pairs (77, ^) for which the pair production upper threshold 
~ 1020 eV. 

shows that in this case, for n = 3 one would get \^^^^\, \ri^^^\ < 10^^^, basically ruling out 
this model. Remarkably, also the case n = 4 would be strongly constrained as in this 
case one could deduce l^'-'^-'l, ^ 10~^. 

This would be the first strong and robust limit on the n = 4, CPT even LV QED, 
which, as we explained in the Introduction, section [1], is also favored from a theoretical 
point of view. Accidentally, as a methodological remark, this result also shows that 
while much attention was focused on the detection of single events at GZK energies for 
constraining LV, full spectral information could be more effective. 

In conclusion, we think that these results fully show the crucial role that the PAO 
experiment will play in coming years as the main experiment for testing fundamental 
symmetries of nature to unprecedented precision levels. 
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